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Abstract. We consider systems of particles coupled with fluids. The parti- 
cles are described by the evolution of their density, and the fluid is described 
by the Navier-Stokes equations. The particles add stress to the fluid and the 
fluid carries and deforms the particles. Because the particles perform rapid 
random motion, we assume that the density of particles is carried by a time 
average of the fluid velocity. The resulting coupled system is shown to have 
smooth solutions at all values of parameters, in two spatial dimensions. 
Key words Nonlinear Fokker-Planck equations, Navier-Stokes equations, 
microscopic variables, Deborah number. 
AMS subject classification 35Q30, 82C3f, 76A05. 

1 Introduction 

We discuss global regularity of solutions of systems of equations describing 
fluids with particle suspensions. The particles are parameterized by indepen- 
dent microscopic variables m that belong to a compact, connected, smooth 
Riemannian manifold M of dimension d. Derivatives with respect to the 
microscopic variables are designated by the subscript g. The particles are 
included in a fluid in IR n , n = 2 obeying the forced Navier-Stokes equations. 

The forces exerted by the particles on the fluid are expressed through 
the divergence of an added stress tensor. The added stress tensor r p (x,t) is 
obtained after averaging out the microscopic variable and the Navier-Stokes 
equation is macroscopic. The microscopic inclusions at time t and macro- 
scopic physical location i G I 2 are described by the density f(x,m,t)dm 
where dm is the Riemannian volume element in M. The density is nonneg- 
ative, / > and 




(1) 



holds for every x, t > 0. 

The added stress tensor is given by an expansion 



oc 




(2) 



fc=i 



where 




(3) 
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T-f(x,t) = / -flf(m 1 ,m 2 )f{x,m 1 ,t)f(x, m 2 ,t)dmidm 2 , (4) 



M J M 



and, in general, 



J ^ i j\m 1 ,...,m k )f(x,m 1 ,t)f(x,m 2 ,t)---f(x,m k ,t)dm ( " 

Mx—xM 

Expansions of this kind for the added stress tensor r p are encountered in the 
polymer literature ([HI)- In it was proved that only two structure coef- 
ficients in the expansion, 7y , 7^ are needed in order to have energetically 
balanced equations, provided certain constitutive relations are imposed. The 
energy balance confers stability to certain time-independent solutions of the 
equations. In this work we are interested only in general existence results, 
and do not need to use special constitutive relations. We will only use the 
fact that the coefficients 7}- are smooth, time independent, x independent, 
/ independent. When infinitely many coefficients are present, we will use a 
finiteness condition assuming that the series 

00 

y^k 3 \\j>^\\ H Pk(Mx-xM) (6) 

fc=l 

converges for a sequence p k > fc+4 2 d+6 . 
From d El ED ) it follows that 

\T P (x,t)\ < cp(x,t) (7) 

holds with a constant that depends only on the coefficients 7^. The spatial 
gradients of r p are of particular importance for regularity. The fact that the 
constitutive coefficients 7^- are smooth functions of the microscopic variables 
allows us to relate the size of the spatial gradients of r p to a rather coarse 
average on M: differentiating (J3J) with respect to x it follows from © that 

\V x r p (x,t)\<cN(x,t) (8) 

(k) 

holds with a constant c that depends only on the smooth coefficients 7^ . 
Here 

N(x,t) = \\RV x f\\ L 2 {M) 
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and 

with s > 1 + 1. The inequalities (0) and (JHJ) are the only information con- 
cerning the relationship between t p and / that we need for regularity results 
in this work. We can use the detailed form (J2J) and the fmiteness condition 
(jBJ) to deduce them, but we could just as well require them instead of (J2J). 

The particles are carried by the fluid, agitated by thermal noise and 
interact among themselves in a mean-field fashion, through potentials that 
depend linearly and nonlocally on the particle density distribution / ([15J). 
We assume that the fluid does not vary much in time during a characteristic 
relaxation time of the particles. Mathematically, this means that the particles 
are carried by a short time average of the fluid velocity. This assumption 
alllows us to prove global existence of smooth solutions and to bound a priori 
the size of the physical space gradients of the stresses. 

The mathematical study of complex fluids is in a developing stage. Most 
results are for models that are macroscopic closures, that is, in which t p 
has its own macroscopic evolution, coupled with the fluid: the microscopic 
variables do not appear at all. Existence theory for viscoelastic Oldroyd 
models is presented in ([IB]); see also ([IB]) for related issues. There are 
few other regularity results concerning complex fluids, including some that 
retain microscopic variables. Among them are ([Zj, ^Oj, [TT] . [14*] . |17j). For 
Smoluchowski equations coupled with fluids, the case in which u is given by 
a time independent linear Stokes equation in n = 3, M = S 2 with r p given 
by a relation © was studied in (jTHj) for the case of a linear Fokker-Planck 
equation and in (0j) for general nonlinear Fokker-Planck equations. 

The proofs in the present work are based on a few key facts. The first 
one is that gradients of bounded by N, and N is controlled linearly 

by the advecting velocity, taking advantage of the a priori boundedness of p 
in L 1 fl L°°. This idea was used in ( 4J) to prove regularity for the system of 
particles coupled with the Stokes system for the fluid. The second significant 
fact concerns the Navier-Stokes system driven by the divergence of bounded 
stresses. We are interested in the size of the time integral of the supremum 
of the norm of the gradient of velocity. This is an important nondimensional 
magnitude that controls the amplification of gradients of passively advected 
scalars. We prove a logarithmic bound for this amplification factor. The 
strategy of proof uses a natural idea introduced in ([3]): time integration is 
performed first in each wave-number shell, to take advantege of the rapid 
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smoothing of small scales due to viscosity. 



2 Estimates for 2-D Navier-Stokes Equations 

Consider a Navier-Stokes system in R 2 : 



where a(x, t) is a symmetric two-by-two matrix that, in this section, will be 
considered to be a given function. We will be interested in estimates when 
a is integrable and bounded by constants that are known a priori and are of 
order one. The physical space gradients of a are possibly large. The aim of 
the bounds is to find the effect that these gradients have on the stretching 
amplification term 



— — h u ■ Vu — Au + Vp = 

at 

V • u = 



V-<r, 
0, 



(9) 




We take smooth, divergence-free initial velocities 



u(0) e L 2 (R 2 ) n W 2 ' r (R 2 ) 



with localized vorticity u = V • u, 



cu(0) E L 2 (R 2 ) nw 1>r (R 2 ) 



with r > 2. We recall the energy estimate 



sup \\u 

t<T 




and the fact that in two dimensions the vorticity obeys 



d t uJ + u ■ V x u — Au = V x ■ div x a. 



(11) 
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Lemma 1 Let r > 2. There exists a constant c r such that 



sup||V«(t)||£ P(<fa) < c r I ||Vt7(t)||| r tZt + c.||u;(0)||| r 



t<T JO 

holds. For r = 2 we have additionally, 



(12) 



T 



\\Au(t)\\ 2 dt < c 2 / \\Va\\ 2 LHdx) + \HO)\\h- (13) 

'0 JO 

Proof. We multiply (fTTJ) by oj r ~ 1 , r > 2 and integrate by parts to obtain 
^Mli^dx) + J u r ~ 2 \V x u\ 2 dx < (r - 1) J \V x a\\V x u\uj r 



d 

rdt 1 



,r-2 



We use a Holder inequality with exponents r, 2, and then with exponents 



2,2: 



< (r-l)\\V x -a\\ L r 
(r- l) 2 



< 



uj r 2 \S7 x u)\ 2 dx 
1 



1/2 



r-2 



Vx • ^IlL-IMI^ 2 + -( / uu r ~ 2 \V x cu\ 2 dx) 



(14) 



which implies that 



sup \\u{t)f Lr <C(r-l) 2 / || V a 
te[o,T] Jo 



■cr(t)||£ r dt + ||u;(0)||i 



(15) 



and thus (|12|) follows. In order to obtain (|13|) we integrate ()14|) in time at 
r = 2. 

We need a logarithmic inequality for Such inequalities were first 

introduced in ([T]). We will write log* (A) = log (2 + A) for A > 0. Note that 
log,(0) > and log, (A//) < log, (A) + log„(/x) holds for A > 0, li > 0. We 
check our inequality 



K*)IU~ < c r \\u(t)\\ L 2 { i + 



\ 



log, 



U t L r 



Wit V 



Ht)\\ L 2 



(16) 



directly from the Biot-Savart law: 

u(x,t) = / (z -1 ) u(x — Z, t)-r—r 

2tt J R 2 \z\ 

where z = A. We pick two numbers < I < L, take a smooth radial function 
> < a(s) < 1, that equals identically one for s < 1 and identically zero 
for s > 2, and write 



,t) = ^- [ (z^uix - z,t)pr, 

Z7T J\z\<l \ Z \ 



(l) / \ 1 f \ i \ ( \z\ \ dz 

uW(x,t) = — [z )u[x-z,t)a 



and 



Clearly 



Z7r J| 2 |>i V L J \ z \ 

u {L) (x, t) = l- [ (^) - *, t) (l - a (M) ) 



U = + U/^ + W(L) 



holds pointwise. It is also clear that 

\u®(x,t)\ < \Wt)\\Lrl^ 

and that 

\u$(x,t)\<P(t)\\ L Aog* ^j. 
We integrate by parts in the term tim, using to = V -1 • u and deduce 

|«(l)0M)| < C— ||ii(t)||i2 

We choose 

' V(*)ll^ 



and 



if, with this choice, it turns out that I < L. If not, then we still take L as 
above, but we take I = L. The inequality (J16|) follows. 
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Let us consider 

= IHI^ow = \ I Mt)\\ P L r {dx) dt\' (17) 



and 

Bl p) = ||V • a\\ LP{0>T;L r {dx)) = { I ||V ■ a(t)\\ p Lr dty . (18) 



T 



p 



Note that (|TU|) implies that 

r T 

swp\\u(t)\\ 2 L2 + / ||Vw(t)||J 2 ^< K . (19) 

with 

K = (Ki 2) ) 2 + \\u \\ 2 L2 (20) 
and note that (fT2j) implies that 

sup ||w(*)|| L r < c r (B^ + \\u(0)\\ Lr ) = n r . (21) 

t<T 

Finally, using (pHIj) and (jTHj) we can write 

jf |K*)||^ < c(l + T) \\\a\\l* {0>T . m) + \H0)\\m}- (22) 

Now we integrate the square of (fTBj) in time, taking the supremum in 
time of the logarithmic part using (|2ip. and bounding the time integral of 
the sqaure of the gradients using (JTHJ). 

Lemma 2 For r > 2 i/iere exists a constant c r such that 

[ \\u{t)\\ 2 L oodt < c r K {log,(fi r ) + log*(fi 2 ) + log.^o} (23) 
./o 



Theorem 1 For r > 2 there exists a constant c r such that 

|Va.Va;U||i,2( 0)T .(£,r)) < 



CrTUVwolUr + c rV ^Q rV /{\og,{Q r ) + log,{Q 2 ) + log,(K )} 

holds. 



(24) 
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Proof. We represent 

VMt) = T i(V^ • <?) - T 2 {uu) + e tA (V^(0)) (25) 
where T\ and T 2 are operators of the form 

hit) ^Th= [ e (i - s)A AM(s)ds 
Jo 

with HI = M{D) homogeneous of degree zero. Such operators are bounded 
in L p (dt; L q (dx)) for 1 < p, q < oo by the maximal regularity of the heat 
equation combined with the boundedness of the operators H in L q spaces 
(see, for example (Q.2J)). The inequality follows then from the bound 

T 

\\u{t)uo{t)\\ 2 L rdt < c r K n 2 r {log„(0.) + log,(0 2 ) + \og,{K )} (26) 

which, in turn, follows from (|2*T|) and (|23|1. 

We mention also an a priori bound for J Q T HitU^oodt for p < 2. This is 
obtained as follows: we write 

uj(t) = e tA uj + 

+ j e (t-s)A | V ± .div x a(s) + d 2 d 1 {ul -u\){s) + {81 - df) Ul u 2 {s)} ds 


(27) 

This follows from a well-known identity 

u-Vu = -d 2 d x {ul - u\) - {8 2 2 - dl) Ul u 2 

An easy calculation verifies this after one writes u\ = — c^Vs u 2 = diip, 
uj = Aip. From (|2Tj) we get, for any 1 < p, q < 00 

IM|i>[(0,t);L9(dx)] < + Cpq\W\\LPl(0,t);Li{dx)] + ^pq 1 1 ^ 1 1 L 2 P [(0,i ) ,L 2 <? (dx)] ( 2 ^) 

We know however from (|19|) that 

||w||L°°[(0,t);i 2 ((ia;)] < Kq 

and 

IK||x2[(0,t);i r (<ix)] < Cr 
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a priori, for any r < oo, with K , C r independent of t G [0, T]. (These 
constants may depend on T because the norm of a in L 2 may depend on T). 
Then, by interpolation 

< C p? (29) 
holds for q > 2 and p < In view of [2E]we get that 

\\^\\LP[(0,t);Lt(dx)] < C pq (30) 

holds for q > 2 and p < Then, taking q > 2 and using a Sobolev 

embedding theorem, we obtain that 



u\\ P L oo {dx) dt < C p (31) 



o 



holds a priori, for any p < 2. 

For the bound for J Q T ||Vw(t)||L°°^ we will use the Littlewood-Paley de- 
composition. Let T>(Q) denote the set of C°° functions compactly supported 
in Q. Let C be the annulus centered at 0, and with radii 1/2 and 2. There 
exist two nonnegative, radial functions x an d f, belonging respectively to 
V(B(0, 1)) and to V{C) so that 

X(0 + S i > ^2-^) = 1, 

and 

\j - k\ > 2 =^ sup(v?(2^')) n sup(<^(2~ fe )) = 0. 

We denote by the Fourier transform on M 2 and let fa, ft., Aj, Sj(j e N) 
be defined by 

ft = F~ l ip and ft = T Xt 
S> = ^- 1 ( X (2-^)^«) = 2 2j J h(Vy)u{x - y)dy. 

Then 
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u = Sou + A?(ii) 

where u G S', the space of tempered distributions, and the equality holds in 
the sense of distributions. 

The well-known Bernstein inequalities (see, for instance j2j) express the 
fact that Aj is localized around the frequency 2K 

Proposition 1 (Bernstein Inequalities) Let a, b G [l,oo], j > 0. There 
exists constants c, independent ofa,b,j such that the following hold: 

WAje^u-WLa < ce-* 22U " 13 \\A jU \\ L », (32) 

IIA.-/IU- <c||/|U-, (33) 
and, in two space dimensions: 

||5,/||^<c(||/|| L2 + V7||V/|| i2 ) (34) 

and 

HA^ulU- < c2^ +2 ^ 1/b - 1/a) \\A jU \\ Lb (35) 
where \a\ is the length of the multiindex a. 

The Littlewood-Paley decomposition is best suited for Besov spaces B* de- 
fined by requiring the sequence 2 SJ || Aj(«)|| L p to belong to £ q and by requiring 
Sq(u) to be in LP . L 2 based Sobolev space norms can be computed in terms 
of the Littlewood-Paley decomposition: 

\\u\\h~\\So(u)\\l. + J2 22SJ \\ A M\l* 

where ~ means equivalence of norms. However, the norm we are interested 
in is the L 1 (0,T;W 1,o °) norm. The C s norms can be computed as 

\\u\\c° ~ ||£q(u)||c« + sup2 JS ||A i (M)|| L oo 

j>0 
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but only if s is not an integer. In order to obtain L°° bounds for the gradient 
we will have to resort to the inequality: 

||Vu||l~ < ||Sb(Vu)|| L <» + Yl ll A i(VM)|| L oc. 

This inequality reflects the embedding B^ 1 C L°°, which is a strict inclusion. 
The advantage of using this sum (the norm in B^ x ) is that we can commute 
time integration and summation, while time integration and supremum do 
not commute in general. 



Theorem 2 Let u be a solution of the 2D Navier-Stokes system (OJ), with 
divergence-free initial data u E iy 1,2 (IR 2 ) n VF 1,r (R 2 ). Let T > and let the 
forces V • a obey 

a E L 1 (0,T;L ao (R 2 )) H L 2 (0, T; L 2 (M 2 )) 

and 

V- a E L^O, T; L r (R 2 )) H L 2 (0, T; L 2 (R 2 )) 

with r > 2. There exists a constant c depending on r such that, for every 
e > 



J Q \\Vu\\ L °°dt < cVT\\u(0)\\ H i + 



+ e 



(36) 



c [k^T + K® log, (^) + (1 + T)Ei log, ( 
holds, where 

Ei= [ \Ht)\\mdt. (37) 



and 

«i = {lHlWi^) + H°)llHi}- ( 38 ) 
Consequently, in view of and Uh}) a&ove 

/ ||V^t)|| L ~^ <c(l + T) 2 inog,(5) (39) 

mt/i K = K + + K% + ||it(0)||#i depending on norms of a and the 

initial velocity, but not on gradients of a, and B = Br + Ri depending on 
norms of the gradients of a. 
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Remark. The bound is in fact for the stronger norm of u in the inhomoge- 
neous space L 1 (0,T; 1 ). 

Proof. We start with the Duhamel formula for the gradient of solutions of 
© 

Vu = e tA Vu + J*e( t - s '> A AU(D)(u(s)®u(s))ds 
- J^e^ A Am(D)(a(s))ds 

The homogeneous operator W(D) is given by 

U(D)(u ® u)ij = Rj(5u + RiR^Rkimuk) (41) 

with Rj = dj(—A)~^ Riesz transforms. The strategy is based on an idea 
of Chemin and Masmoudi ([3 ) to take the time integral first, for each wave 
number shell. They did not use information about derivatives of a, and 
therefore obtained only bounds for sup g>1 J Q T || AqVu\\L°°ds. We will use the 
gradients of a to bound the high frequencies and will sum in q in order to 
estimate || Vw||i,i((o,t);i 00 )- Also, we use somewhat different estimates than 
them for the individual shell contributions, but like them, we take advantage 
of a time integration at each shell. We treat separately the contributions 
coming from a and those coming from u®u: 

xj u = F + U + e tA Vu 

where 

F(t) = - [ e {t - s)A AM(D)(a(s))ds. (42) 
Jo 

and 

U(t)= [ e {t - s)A AU(D)((u®u)(s))ds. (43) 
Jo 

Clearly 

/ \\SoF{t)\\ L °°dt < cK^T (44) 
Jo 

is true using for instance \\S (F(t))\\ L oo idx) < - A)S (F(t))\\ L 2 {dx) . We 
take q > and apply A q : 

||A 9J F(t)|| L « < /V^^V^IHs)!!^ (45) 
Jo 
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(40) 



Integrating on [0, T] and changing order of integration we obtain 



A q F(t)\\ L o,dt< / Ks)|M/ e 



-(t-s)2 2 <9- 1 ) 



T 



2 2q dt)ds 



<c \\a(s)\\ L oods. (46) 
Jo 

We bound the same quantity differently, with large q in mind, and use 
Vtr G L r with r > 2: 



lA^IU- <c / e -2><«-i)(t-.) 2 ,(i+2/r)|| V . allLrds 



(47) 



We integrate on [0, T] and change the order of integration as above, to obtain 



T ||A,F(*)|| £ <»di < f || V ■ a{s)\\ Lr { f e- 22(q ~ 1)(t - s) 2^ 1+2 ^dt)ds 

JO Js 

T 



< 



C 



2«(l-2/r) 



\V ■ a(s)\\ L rds (48) 



o 



Using (jUJ), (jinj) to estimate the small wave numbers in F and (pfSj) to 

estimate the high ones, we obtain 



£\\F(t)\\ L ~dt< 

J T [\\S F(t)\\ L ~ + J2o< q <m \WF{t)\\w + E 9 >m \\\F(t)\\L~] dt (49) 
< cK^T + cM J T H^IUoods + 2M(1/ c 2 _ 1/r) J T || V • ff(s)|| L rds. 

Then, choosing M 



c Io ' H V " a ( s )\Wds 



M = c r log, 



(50) 



we obtain 



\\F{t)\\ L ^dt<cK^T + cK^ log 



B 



(i)' 



+ e 



(51) 



with AT^ defined in (J]3j) and B} ij defined in JEJ. 



(i> 
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Remark. We do not need to integrate the term F in time, if a is bounded: 
we can obtain a pointwise logarithmic bound for F(t) in terms of Br with 
p > Indeed, from we have 

||A g F(t)|| L oc < c\\<j\\ L ™( dt dx) (52) 

and from (JTTj) we obtain 

\\A q F(t)\\ L oo < c2- 2q ^-- P ) \\W ■ a\\ LP(0tT . Lr) (53) 

Summing (|52|) from q = to q = M, summing (|53|) from q = M to infinity 
and choosing M appropriately, we obtain 

/ b (p) \ 

sup \\F(t) \\ L oo <C^T + cK 00t00 log, — — (54) 

t<T \ ^00,00 / 

with 

K<x>,ce — || cr ||i: oo (0,T , ;i 00 ) ■ (55) 

This bound can be used to reprove the global existence of the Stokes system 
coupled with nonlinear Fokker Planck equations. 
We split the nonlinear term 

U(t) = S 2 (U)(t)+V(t) (56) 
S 2 (U) = S (U) + A 1 (U) + A 2 (U), (57) 
V(t) = J2\(U) (58) 

q>3 



with 
and 

Clearly 



/ \\S2(U)\\L™dt < cTE 1 , (59) 
Jo 

where E\ is given in (|57j) . For the nonlinear term V we use Bony's de- 
composition (see for instance j2]) into commensurate and incommensurate 
frequencies: 

V(t) = C(t) + I(t) (60) 
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with 



C(t) = J2 e ( '- s)A Ai(D)A, \{u{s)) ® A p ,(u{s)) \ ds (61) 

<7>3 J ° \|P-P'|<2 / 



and 



J (*) = X)/ e (t ~ s)A AH(L>)A 9 ( ^ A p (u) <g> Ap/(u) I cfe (62) 

g>3 ^° \b-p'|>3 / 

In the decomposition above we made the convention that the indices p,p' 
run from —1 to infinity and when p = —1 instead of A p we have Sq, and of 
course, the same thing for p'. 

We start by treating the term C(t). Because the range of q is q > 3 it 
follows that the range of p,p' is p > 1, p' > 1. Then we estimate inside the 
integral 

\\A q {{A p u{s)) ® (Ap/(u(s)))||x,oo < c2^\\A q ((A p u(s)) ® (A p ,(u(s)))\\ L i 

< c2 (2 9 -2 P )|| VApM(s) || L2 || VAp/M(s) || i2 

Clearly at least one of p,p', say p, satisfies p > q—2. Now, because \p— p'\ < 2, 
it follows that p' = p + j, j 6 [—2,2] and therefore, changing the order of 
summation in (J61j) 

\\C(t)\\ L ~ < E, 2 =- 2 E P °°=i c /o l|A p Vn( S )|| L2 ||A p+J Vn( S )|| L2 
x |^3 2 e- 22( ^ 1) ^ s )2^2 2 ^)} ds 

Integrating in t and changing order of integration we obtain 

\\C(t)\\ L ~dt < 

cE J U 2 Er=Jo T ll A P ( VM ( s ))IUH|A P+J <v M ( s ))|| L2 ES2 2(9 - p) ^ 

< cYS^IoT^i \\A p (Vu(s))\\ L 4A p+j (Vu(s))\\ L2 ds. 



We have obtained thus 

•T 



[ \\C(t)\\ L oodt < cE v (63) 
Jo 
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We turn now to the term I(t) of (|62j) . This term is made up of two sums, 



I(t)=h(t)+I 2 (t) 
h(t) = E g > 3 Si e (< - s)A Ai(D)A 9 (E P > p/+ s \«s)) ® A p ,(u(s))) ds 
h{t) = E g>3 Jo e (< - s)A Ai(D)A 9 (£ p , >p+3 ® M«00)) 

(64) 

We will treat ii because the treatment of /2 is the same, mutatis mutandis. 
Because p > p' + 3 it follows that, in order to have a nonzero contribution at 
q, the index p must belong to [q — 2, q + 2], i.e., p = q + j with j G [—2, 2]. 
Then we can write 

Jl W =EE / e ( '- s)A AH( J D)A g (J,( S ))^ (65) 

with 

J 9 (s) = A q+j (u(s)) ® S, +i _ 3 (it(s)). (66) 
For q < M we estimate 

||A 9 (J ? ( S ))|U~ < c||^ +i _ 3 (K( S ))|Uoo||A g+J ( w ( s ))|| J 



(67) 



< c [\\u(s)\\ L 2 + VM + 2\\\7u(s)\\ L2 ] \\A q+j (u{s))\\ La o 
where we used f|34)l : 

||5, +i («(a))|Uoo < c (||n( S )|| L2 + v/^+JllV^s)!^) . (68) 

Using Bernstein's inequality || A q+ ju(s)\\i,°° < c||A g+J -Vw(s)||i2 we obtain 



l|A g (J 9 (s))|| L o= < c (j|u(s)|| L2 + y/M + 2\\Vu(3)\\v) \\A q+j V(u(s))\\ L , 

(69) 

For q > M we estimate 

l|A 9 (J 9 00)|U~ <c(|K5)|| i2 + ||A U ( S )|| i2 )2^||A g+i AK S ))|U 2 (70) 
We write 



/ \\h(t)\\ L oodt<A + B 
Jo 
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with 



2 JSd 'j 1 / T \ 

A = cJ2Y,[ \\\(J q (s))h~ ( [ 2^e-^ 2iq - 1) dt)ds 

j=-2 g=3 J ° ^ Js ' 



and 



B = cJ2H I W A Ms))h~ ( I 2^e-^ 2(q - 1] dt) ds 

j=-2q=M J ° ^ Js ' 

We use dHHD for A: 

A<ctf (\\u(s)\\ L2 + VM + 2\\Vu(s)\\ L ,) (e^ 3 ||A g+i V«(5)|U 2 ) 



< c / Q T (||w(s)|| x2 + > /M + 2||Vu(s)|| ra ) (VM||Vu( S )|| L2 J ds. 
We have therefore 

A < cMEi (71) 

We use (J70D for 5: 

5 < c / T (||«(5)|U 2 + ||A U ( S )|| i2 ) (2- q \\A q+J A(u(s))\\ L2 ) ds 

JO „^»r 



q>M 



and therefore 



where 



B < c2- M f \\u(s)\\ 2 H2 ds = c2~ M E 2 
Jo 

E 2 = f \Ht)\\ 2 H2 dt. (72) 
J o 

In view of (J22J) 

B < c2- M (l +D {lMli 2(0 , W + ll«(0)|||i} = c2- M (l+T)R 1 (73) 
For any e > 0, we choose 

t , , /c(l + T)iV 
M = log, ' 
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with Ri given above and we obtain from (|7Tj) and (J75j) 

||J(t)||x-<ft < cfrlog. + p Rl ^j +e (74) 

The sum of the inequalities (J5TJ) . (|o^jl. (jUj) and a straighforward esti- 

mate for the linear term carrying the initial data give the inequality (}36|) of 
the theorem. 



3 Coupled Nonlinear Fokker-Planck and 
Navier-Stokes Systems in 2D 

We consider now the coupling between fluid and particles. The evolution of 
the density / is governed by a nonlinear Fokker-Planck equation 

dtf + v- V x f + div g (Gf) = -AJ (75) 

r 

The coefficient r > is the time scale associated with the particles. The 
microscopic variables m are non-dimensional. The tensor G is made of two 
parts, 

G = W g U + W. (76) 

The (0,1) tensor field W is obtained from the macroscopic gradient of velocity 
in a linear smooth fashion, given locally as 

/ - dv \ 

W(x,m,t) = (W a (x,m,t)) a=1 _ d = ^(m)-^(x,t) (77) 

w=i aXj J a=1 . 



The smooth coefficients c^(m) do not depend on the solution, time or x 

(*: 

Yij 



and, like the coefficients 7^, they are a constitutive part of the model. The 



potential U is given by 

U(x,m,t) = -(JCf)(x,m,t) (78) 

where b is a nondimensional measure of the intensity of the inter-particles 
interaction. The nonlocal microscopic interaction potential 

(£/) (x, m,t) = J K(m, q)f(x, q, t)dq (79) 

Al 
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is given by an integral operator with kernel K(m, q) which is a smooth, time 
independent, x independent, symmetric function K : M x M —> M. Q15J). 
The Navier-Stokes equations are 

d t v + v ■ V x v + V x p = vA x v + Va • t p (80) 

with V x -v = 0. The added stresses are given by the relations (J2J), with (jSJEJ)- 
The added stresses are proportional to kT (where k is Botzmann's constant 
and T is temperature) and have units of energy per unit mass. The density 
of the fluid is normalized to one. The particles are advected by 



1 



T 



v(x,t) = — v(x,s)ds (81) 



(*-t)+ 



We rescale the Navier-Stokes equations using the length scale A = v/^ and 
the time scale \ 2 v~ l = Z, where 5 is the Deborah number, the ratio between 
the relaxation time scale of the particles and the macroscopic (observation) 
advective time. We set 

, . .A fx t5 
v[x,t) = b-u -, — 
r \ A t 



and we arrive at Q with 



a = T -r r . (82) 



The Fokker-Planck equation becomes 



{dt + u- V x ) f + div fl (G/) = ^A g f (83) 



with 

ft 



l r 

u(x,t) = - / u(x,s)ds, (84) 
6 J(t-s)+ 

and (|7B|) with 



n du \ 

W(x,m,t) = (W a (x,m,t)) a=1 _ d = ^ ^( m )^(x,t) (85) 

and 

U(x,m,t) = -{JCf)(x,m,t) (86) 
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The forces applied by the particles are obtained after / is integrated 
along with smooth coefficients 7^ on M in order to produce a. Therefore, 
only very weak regularity of / with respect to the microscopic variables m 
is sufficient to control a. In order to take advantage of this, we consider the 
L 2 (M) selfadjoint pseudodifferential operator 

R=(-A g + iy% (87) 

with s > I + 1. We will use the following properties of R: 

[R,V x ] = 0, (88) 



RV g : L\M) -> L 2 (M) is bounded, (89) 

RV g : L 2 (M) -> L°°(M) is bounded, (90) 

[V g c, R- 1 } : H S (M) -> L 2 (M) is bounded, (91) 
for any smooth function c : M — ► R, and 

i? : L 2 (M) -> i? s (M) is bounded. (92) 

We differentiate (}83|) with respect to 2, apply R, multiply by RV x f and 
integrate on M. Let us denote by 

N(x,t) 2 = J \RV x f(x,m,t)\ 2 dm (93) 

M 

the square of the L 2 norm of RV x f on M. The following lemma was proved 
in (jH): 

Lemma 3 Let u(x, t) be a smooth, divergence-free function and let f solve 
A83\). There exists an absolute constant c > (depending only on dimensions 
of space, the coefficients and M, but not on u, f , 5) so that 

(d t + u- V x ) N < c{\V x u\ + -)N + c\ V x V x u\ (94) 



holds pointwise in (x,t). 
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The proof is given below in the Appendix for completeness. It works inde- 
pendently of the dimension n of the variables x. The equation obeyed by 

P = Im f dm is 

{d t + u- V x )p = 0. (95) 
We will take initial densities that obey 

< p(x,0) < 1. 

Therefore 

0<p(x,t)<l (96) 

continues to be true and, in view of (J7J) and the fact that u is divergence free, 
it follows from (J55J) that 

Mt)\\ L r <cv (97) 

holds if we assume (as we do) that p(x, 0) G L : (R 2 ). The inequalities (0) 
and ((HI) use only the smoothness of the coefficients 7, the relations (0 0)) 
the condition (JBJ) and (0), and therefore they hold throughout the evolution. 
We know thus that a is bounded in L°° and that K^ p) (in JQ) are bounded 
a priori. (We also know that J Q T ||w(£)||£,°od£ is bounded a priori (see (|3*T|)). 
This implies that the support of p, if initially compact, would expand only a 
finite amount in finite time. We do not use this for the proof, but obviously 
this is a physically important a priori quantitative information.) 

Theorem 3 Consider the coupled Fokker-Planck and Navier-Stokes system 
|75| ), / fffflj) wift arbitrary parameters v,T,b > 0. Assume that the initial 
velocity v(0) is divergence-free and smooth, v(0) G W 2 ' r (R 2 ) n L 2 (M 2 ) 
r > 2. Assume that the initial distribution of particles f(x,m,0) is non- 
negative, smooth, in the sense that 

N(x,0) = \\V x f(x,;0)\\ H - s[M) G L r (R 2 ) fl L 2 (R 2 ) 

for some s G 1R, and localized, in the sense that 

p(x, 0) = / f(x,m,0)dm 
Jm 

obeys < p(x, 0) < 1 and p(-,0) G L 1 (M 2 ). T/ien the solution of the system 
d7<5[ f^) exists for all time and is smooth. In particular, the norms of 

v G L\dt- W 1 ' 00 ^ 2 )) n L 2 {dt; W 2 ' r {E. 2 )), 
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/ G L°°(dt; W^°°{dx] H- S (M))). 

can be bounded a priori in terms of the initial data, for arbitrary large finite 
intervals of time [0, T] . 

Proof. Let T > be given. The sort time existence of solutions and the 
uniqueness of the solutions can be obtained following classical methods of 
proof. Let 

n(t) = \\V x a\\ 2 Lr + \\a\\ 2 Hl , (98) 



B{t) = [ \\V x a\\ 2 L rdt+ [ \\a\\ 2 H idt, 
Jo Jo 



(99) 



and 

Using (J9"2j) we have 



0(t) = ||v«(t)|U«, (ioo) 



7 (£) = sup g{s), (101) 

0<s<t 



G{t) = \\V x V x u(t)\\lr + \\u(t)\\h- (102) 



^ < cg(t)n(t) + cG(t) (103) 
dt 



Integrating on (0, t) 

n{t) < n(0) + cj(t)B(t) + c [ G(s)ds (104) 

Jo 

In view of (136)1 we have 

crf(t) <G){l + log*(£K*))} (105) 



with Co a constant depending on T and the initial data. In view of ()22|) and 
P4"|) we know 

f G{s)ds<C l {l + B{t)\og,B{t)} (106) 
Jo 

with C\ a constant. Note that 

n(t) = jB{t) (107) 
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Thus, from (|TU4l . (ITUo]) and (ITUfiJ) we deduce 

jB{t) <C 2 {1 + B{t)\og,B{t)} (108) 

holds with C 2 a constant that depends on T and the values || Vu;(0)||i/-, 
||o;(0)|| i 2, ||tt(0)||i2 and n(0). This produces a pointwise-in-time a priori 
finite bound for B(t) on the interval [0, T], and retracing our steps, via (jl06|) 
and ()104J) . on n(t). Once the forces in the two-dimensional Navier-Stokes 
equations are known to be thus bounded, it follows (from (|3Tj) -but also much 
easier, from energy estimates) that the solution of the Navier-Stokes equation 
is smooths as stated. 



4 Appendix: Proof of Lemma 3 

The evolution equation of N is 

- (d t + u ■ V x ) N 2 = -D + / + // + III + IV (109) 



with 



and 



D = e J \V g RV x f\ 2 dm (110) 



M 



du 3 f ( R df_)( R dL ]dm (m) 



dxk J V & X j J \ ® X k 

M 



II = -E( V ^) / (Rdiv g (<%f))(V x Rf)dm, (112) 

III = - E ^ / (Rd™ 9 (c%VJ))(RV x f)dm, (113) 



a=1 ~ M 



IV = -^J Rtivg(V t {fV g (JCf)})RV x fdm. (114) 



M 



Now we start estimating these terms. We will use repeatedly and 
(|96|). In view of the fact that D > 0, we may discard this term. Clearly 

|/| < c\W x u\N 2 . (115) 
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In order to bound II we use (|83j) to bound 

\\R^ 9 {4f)\W m <c\\f\W {M ) = c 

so that we have 

\II\ < c\V x V x u\N. (116) 

In order to bound III we need to use the commutator carefully. We start by 
writing 

IMW.^Y,./) = Rdiv g (^R- l R\/J) = 
dw g (c^RV x f) + [Rdxvg&Rr 1 ] RVJ. 
The second term obeys 

|| [iMiv^iT 1 ] RV x f\\ LHM) <cN 

because, in view of (|91|) and ()92)) one has that 

[Rdiv g (%, R- 1 ] : L 2 (M) -> L 2 (M) is bounded. 

The first term needs to be integrated against RV x f and integration by parts 
gives 

J (div g (c%RVJ))RVJdm = l -J (dw g ^)\RV x f\ 2 dm. 

M M 

We obtain thus 

\III\<c\V x u\N 2 (117) 
The term IV is split in two terms, IV = A + B 

A = -- 5 j Rdiv g ({(V x f)V g (}Cf)})RV x fdm (118) 



M 



and 



B = -- 5 j Rdw g ({fV g (JCV x f)})RV x fdm. (119) 



M 



The (0,1) tensor <&(x, m, t) = (V g JCf)(x, m, t) is smooth in m for fixed x, t 
and 

\\${x, •,t)||w i .°°(Af) < c s 
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holds for any s, with c s depending only on the kernel K. We write the term 
A 



A = -I J Mv s ({(VJ)$})i?VJd 



m 



M 

_ /' T 



b J R- l (RV x f){$ ■ V g R 2 V x f))dm 

M 

= -±1 div 9 {$} \RVJ\ 2 dm+IJ (RVJ) [R~\ $V 9 ] R(RV x f)dm. 

M M 

In view of (JEJ), the operator 

[iT 1 , $V g ] i2 : L 2 (M) -> L 2 (M) 
is bounded with norm bounded by an a priori constant. It follows that 

\A\ < ^N 2 (x,t) 
o 

holds. The term B is easier to bound, because 

()CS/ X f)(x, m, t) — J R~ 1 K(m,n)RV x f(x,n,t)dn 

M 

and thus 

\\{VglCV x f){x,;t)\\L°°(M) < cN{x,t). 

Using @2| it follows that 

\B\ < ^N 2 (x,t) 
o 

and consequently 

\IV\ <^-N 2 (x,t). (120) 
o 

Putting together the inequalities (jll5j) . ()116|) . (|117|l and (|120jl we finished 
the proof of the lemma. 
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